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It has been shown that the Schwinger-Dyson equations for non-Hermitian theories 
implicitly include the Hilbert-space metric. Approximate Green functions for such 
theories may thus be obtained, without having to evaluate the metric explicitly, by 
truncation of the equations. Such a calculation has recently been carried out for var- 
ious PT-symmetric theories, in both quantum mechanics and quantum field theory, 
including the wrong-sign quartic oscillator. For this particular theory the metric is 
known in closed form, making possible an independent check of these approximate 
results. We do so by numerically evaluating the ground-state wave-function for the 
' equivalent Hermitian Hamiltonian and using this wave-function, in conjunction with 

• the metric operator, to calculate the one- and two-point Green functions. We find 

D . that the Green functions evaluated by lowest-order truncation of the Schwinger- 

Dyson equations are already accurate at the (6-8)% level. This provides a strong 
justification for the method and a motivation for its extension to higher order and 
to higher dimensions, where the calculation of the metric is extremely difficult. 

PACS numbers: 02.30.Cj, 02.30.Mv, 03.65. Ca, 11.30.Er 

I. INTRODUCTION 

The wrong-sign quartic oscillator, as we shall explain below, is a particularly interesting 
j — . ■ example of a quantum theory whose Hamiltonian is not Hermitian, but is instead PT- 
symmetric, and as a consequence of that unbroken PT symmetry possesses a completely real 
energy spectrum. The recent surge of interest in this type of theory began with the pioneering 
paper of Bender and Boettcher[l|, in which they showed, by numerical and asymptotic 
O 1 analysis, that the entire class of Hamiltonians 
O ' 



H = p — {ix) M (1) 



had that property for M > 2. Apart from the trivial case M — 2, the simplest example 
is for M — 3, where H = p 2 + ix 3 , for which the Schrodinger differential equation can be 
analyzed on the real x axis. However, for M > 4 the eigenvalue equation must instead be 
formulated in the complex plane, inside a Stokes wedge, in order that the wave-function 
vanish exponentially at infinity. In particular, for M = 4 we obtain the wrong-sign quartic, 
H = p 2 — x 4 , which is unbounded below on the real axis, but can be understood as a 
legitimate Hamiltonian when formulated on an appropriate contour in the complex x plane. 

Since that initial paper, there has been intensive investigation into the properties of 
Hamiltonians of this kind, whose progress can be followed in the reviews by Bender j2] 
and Mostafazadeh[3|]. We restrict ourselves here to those features that form the essential 
background to the present calculation. 

For a viable framework of quantum mechanics one needs not only a real spectrum but 
also a probabilistic interpretation. In standard quantum mechanics that is provided by 
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matrix elements of the type J ip*A(f), orthogonality of eigenfunctions J dx ip\ip2 = 0, and 
the probability density In PT-symmetric quantum mechanics it is found instead that 
orthogonality of eigenfunctions takes the nonlocal form j dx ip\(—x)ip2{%) — 0, that is, 
J dx {iPi)pt4>2 — 0. This is a problem because the metric involved in the corresponding nor- 
malization integral J dx (ip) pt4> is n °f positive-definite, and we do not in the first instance 
have a proper probabilistic interpretation. However, it was subsequently found [4] that an- 
other metric could be constructed, with the help of a grading operator C, which preserved 
orthogonality and gave a positive normalization integral j dx (ip)cPT'^- In contrast to stan- 
dard quantum mechanics, this metric is not universal, but is dynamically determined by 
the particular Hamiltonian in question. The calculation of this metric is extremely difficult, 
and in most cases can only be calculated approximately, either through a set of algebraic 
relations jBf, or through the use of Moyal brackets 0. 

A more general framework, of which CPT-symmetry is a special case, was developed by 
Mostafazadeh[3]. A Hamiltonian H is said to be quasi- Hermit ian if it can be related to a 
Hermitian Hamiltonian h by a similarity transformation: 

H = p-'hp, (2) 

where p is a positive-definite Hermitian operator. From this we immediately see that 

H ] = 7]Ht]-\ (3) 

where rj = p 2 . The connection with the CPT formulation is that rj can be identified as 
when CP is written[5] in the exponential form CP = e®. Accordingly p can be written as 
p = e~z Q . 

Matrix elements of operators should always be of the form {ipi\rjA\ip2) , i.e. calculated 
with inclusion of the metric. This result is easily obtained by a similarity transformation 
from the corresponding matrix element in the Hermitian theory. Thus, if are the 

corresponding states and a the corresponding operator, 

if>i = e^ Q (fi, A = e* Q ae~^ Q . (4) 

Hence B 

(fi\a\<p 2 ) = (fale-frie-frAefye-frlfo) = {^\e- Q A\^ 2 ) . (5) 
In view of Eq. (jlj), we can also write 

(Vi|e-°A|^) = ((pxle-^Ae^M, (6) 

a relation we shall need below. 

When one is calculating Green functions using path integrals, the Schwinger- Dyson equa- 
tions or Feynman diagrams, where the metric makes no explicit appearance, it is not im- 
mediately clear that the quantities so calculated correspond to matrix elements evaluated 
with the inclusion of the metric. This issue was addressed in Ref. @, with the conclusion 
that this is indeed the case, the essential point being that only when the metric is correctly 
included do the Heisenberg equations of motion and the canonical commutation relations 
take their standard form. This leaves the way open to using non-perturbative truncations 
of the Schwinger- Dyson equations to obtain approximate values for the Green functions of a 
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quasi-Hermitian theory without the need to evaluate the metric, which, as mentioned above, 
is generally extremely difficult to calculate. 

Such a non-perturbative calculation has recently been carried out by Bender [Kj], for 
various quasi-Hermitian quantum theories, and has the potential to be extended to quantum 
field theories, where the calculation of the metric is even more problematic. The method 
involves truncation of the Schwinger-Dyson by setting to zero the connected Green functions 
beyond a certain order, and has been applied, among other theories, to both the standard 
anharmonic oscillator and the wrong-sign anharmonic oscillator with or without a quadratic 
term. This last theory is special among quasi-Hermitian theories in that it is essentially 
the only non-trivial theory where the metric is known exactly [ill]. Hence in this particular 
case we have the possibility of calculating the Green functions numerically, by utilizing the 
similarity transformation to the equivalent Hermitian theory, thus providing a check of the 
accuracy of the Schwinger-Dyson calculation. 

The paper is organized as follows. In the next section we summarize, for completeness, 



the method of Ref. [10j and the results obtained there, specifically the renormalized mass 
and the first two connected Green functions G\ and G2 for the wrong-sign quartic oscillator 
without a quadratic term. In this case we have scaling, and the result is universal. One 
can also easily include a quadratic term, but then one has a one-parameter set of solutions 
depending on the quadratic coefficient. In the following section we show why previous 
calculations of the Green functions, which predated the realization that matrix elements 
should be calculated with the inclusion of the metric, are in fact correct. We then describe 
how to calculate the Green functions using the Schrodinger wave-functions of the equivalent 
Hermitian Hamiltonian h, finding that, even in the simplest truncation, the Schwinger- 
Dyson equations give quite reasonable values. We conclude with a brief discussion of the 
significance of these results. 



II. TRUNCATED SCHWINGER-DYSON EQUATIONS 



In this section we summarize the results obtained in Ref. [10|, which to this order essen- 
tially amount to mean field theory. The starting point is the expression for the vacuum- 



generating functional as a Euclidean functional integral [12 



Z[J] = J [dip] exp 



(lr ( 2<y) 2 + - - J<p 



(7) 



The connected one- and two-point functions, in the presence of J, are respectively defined 



as 



G ' w = ik) laZlJ] (8) 

G,{x,y) = JjLfGiM- (9) 

The Schwinger-Dyson equations are derived from the Heisenberg equation of motion for (p: 

(-dV(x) -m 2 )<p- gtp s (x) = J(x). (10) 
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This becomes an equation for Green functions when we take its expectation value in the 
vacuum, still with J ^ 0: 

( -fl»- m » )GlW - g «!VWW>£ = J (,). (ii) 

Here the ((...)) notation means that the expectation value is to be taken with the appro- 
priate metric, but the precise form of that metric plays no part in the equations that follow. 
Z itself is given by Z[J] = ((0\0))j. 

Taking successive functional derivatives 5/5 J(x) of the identity ZG\{x) = ((0\<p(x)\0)) j 
of Eq. ([8]) we obtain 

ZG\{x) + 3ZG 1 (x)G 2 (x,x) + ZG 3 (x,x,x) = ((0|^ 3 (x)|0) j, (12) 

which allows us to substitute for ((0\(p 3 (x)\0) j in Eq. ffTTj) . to obtain 

(-d 2 - m 2 )Gi(x) - g [Gl(x) + 3G 1 (x)G 2 (x, x) + G 3 (x, x, x)] = J(x). (13) 

Differentiating once more, 5/5 J(y), gives the equation 

(_ 5 2 _ m Z) G2 ( x) _ g [3Gl(x)G 2 (x, y) + 3G 2 (x, y)G 2 (x, x) 

+3Gi(x)Gs(x, x, y) + G±(x, x, x, x)] —5(x — y). (14) 

We now set J = 0, in which case G\(x) becomes a constant independent of x, and G 2 (x,y) 
becomes a function of x — y, so in a change of notation we write it as a function of that 
single argument: G 2 (x — y). Finally we truncate Eqs. (fl3"]) and fTT4l) by neglecting G 3 and 
G 4 , to obtain 

m 2 G x -g[G\ + 3G?iG 2 (0)] = 0, 

(15) 

[-d 2 + m 2 -3g (G\ + G 2 (0))] G 2 (x - y) = 5{x - y). 



Assuming 13| that G\ ^ 0, the first equation gives 



m 2 - g (Gj + 3G 2 (0)) = 0, (16) 



while the second equation gives 



P -M\x\ 



where 

M 2 = m 2 - 3g(G\ + G 2 (0)), (18) 
so that G 2 (0) = 1/(2M). Hence we have 

M 2 = ~2gG\ = -2m 2 + g (19) 
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Let us restrict ourselves to the case m = 0, in which case 14 



M = (3^) 1/3 , 

iG 1 = M (2Q) 
(2c/) : 



G 2 (0) 



1 



2M 



Since in the massless case the results scale with various powers of g, we may as well take 
some specific value. For ease of comparison with the calculations of the following section we 
choose g = 2, in which case 

M = 6 1/3 = 1.81712, 

id = ~M= 0.90856, (21) 
2 

G.(0) = ± = 0.27516. 

Note that for <? = —2, the conventional anharmonic oscillator, the corresponding results are 
M = 3 1 / 3 , G x = 0, G 2 (0) = 1/(2M). 



III. ALTERNATIVE METHODS FOR CALCULATING GREEN FUNCTIONS 

In fact these Green functions have been estimated already, in a paper that preceded the 
introduction of the rj metric [3], in two different ways. 

The first was to calculate ((z n )) in terms of the Schrodinger wave-function as 

t> (' dzibo{z) 2 z n t x 

where we write z instead of x to emphasize that the theory must be formulated on an 
appropriate contour C in the complex z plane. This is the first difference from a conventional 
Hermitian theory. The second is that the expectation value is taken with ipo(z) 2 rather than 
| , 0o(^)| 2 - This can be motivated by the fact that the former, as an analytic function, can be 
analytically continued, whereas the latter cannot, and moreover it is line with perturbation 
theory, which involves ipQ alone, not its complex conjugate. 

However, in the light of what we have subsequently learnt about the metric of such 
theories, it is natural to ask whether Eq. f l22|) does indeed correspond to the expectation 
value ((z n )), which explicitly involves the metric. The answer is in the affirmative because 
of the following identity. 

Recall that from Eq. (SJ) V'o — e3< "Vo ; where <£>o ; the ground-state wave-function of a 
Hermitian Hamiltonian, can be chosen to be real. Thus 

= e" Q ^ , (23) 

since Q(x, y) (see Eq. 02 8|) ) is pure imaginary (and antisymmetric). Thus, in Eq. ( I2"2"j) we can 
take the first factor of ip as ipo = e Q ipQ. Then, under integration by parts, or equivalently 
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using the antisymmetry of Q, we obtain 

dzip z n ip = [ dzip* e- Q z n i) , (24) 



so that Eq. (122]) indeed gives the correct formula for ((z n )), taking into account the metric. 



The second method used in Ref. [151 ]. which can be extended to higher dimensions, was to 
apply the strong-coupling approximation to the functional integral. Here again the metric 
does not enter explicitly, but nonetheless, as explained in Ref. [9(, the metric is implicitly 
involved, and the Green functions so calculated are of the type ((z n )). 



IV. GREEN FUNCTIONS FROM THE EQUIVALENT HERMITIAN 

HAMILTONIAN 

We start with 

H z =p 2 z -Xz\ (25) 

where, as in the previous section, we write z instead of x to emphasize that the theory must 
be formulated in the complex z plane. The connection with the Hamiltonian of Section II 
is that H z = 2Hsd with A = \g = 1. 

There is one particular contour, parametrized by the real variable x, which makes the 
subsequent calculations practicable, namely fill] 

z = -2y/(l + ix). (26) 
This gives an intermediate non-Hermitian Hamiltonian 

H x = -{(1 + ix),p 2 x } - -p x - a(l + ix)\ (27) 
where a = 16A. For this Hamiltonian Q can be found explicitly as 

3 

Q = -^ + 2p, (28) 
and the equivalent Hermitian Hamiltonian h as 

ii 1 

h = e^ Q He^ Q = ^---p + ax 2 . (29) 
4a 2 

Since h is quartic in p but quadratic in x it is natural to work with the Fourier-transformed 
wave-function (pip). The Schrodinger equation for (pip) is 

{- a W-\ p+ ^) m = E ^ p) - m 



Scaling p to p = yy/a we get 



( d 2 1 y 4 \ 

\dy^ ~ 2 VV ^ + a jj¥> = E <P- ( 31 ) 
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Recall that we have chosen A = 1, so that a = 16. 

Before going on to the Green functions let us first check M, which in terms of quantum- 
mechanical energy levels is the gap between the first excited state and the ground state, 
divided by two because H z = 2Hsd, namely M = \{E\ — Eq). The numerical eigenvalues 
obtained from Eq. fl3Tl) are 

E = 1.4771497535 

(32) 

E x = 6.0033860834 , 



in agreement with Bender and BoettcherpJ, who obtained these results by numerically 
solving the eigenvalue equation along a complex contour. Eq. (132]) gives M = 2.26312, 
showing that the estimate 1.81712 from the lowest-order truncation of the Schwinger- Dyson 
equations is not particularly accurate. 

Turning our attention now to the Green functions, we use the result (cf. Eq. <^) that 
the expectation value of any operator O is obtained in terms of the ground-state WF <po(p) 
of h as 

(0) = {Up)\^ Q Oe^\Up)) ■ (33) 

So in particular 

(z) = -2^ (p)|e-^(l + 2x)M Q |^oG»)> , (34) 

where x = id /dp. 

It is in fact easier to calculate 

(z 2 ) = -4(£ (p)\e-&(l + ix)e^\<p (p)) , (35) 

whose S-D estimate is 

(z 2 )sd = G 2 (0) + G\ = -0.5503 . (36) 
We now check this from Eq. (1351) . Note that 



ixe^ Q o (p) = - (e^o(p)) 



p \ d 



2a) dp 



MP) (37) 



So 



(z 2 ) = -^(Up)\p 2 \Up)), (38) 

since (p${p) can be taken real, so that the derivative term vanishes on integration. In terms 
of y, this is 

(z 2 ) = -2(<p \y 2 \£ ). (39) 
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The numerical value, found by evaluating this expectation value with the numerical ground- 
state wave-function <p(y) is (ip \y 2 \(p ) = 0.2585. Hence 



-0.517, 



(40) 



which is actually quite close to the Schwinger- Dyson estimate -0.5503 of Eq. ( 136 



Returning now to the calculation of (z) from f[3"4l) . we deal with (1 + ix) 2 by writing it 



as 1 



(1 + ix) 2 = (1 + ix)(l + ix) 2 

. . du 

[l + tx) 



-u 2 (1+ix) 



d_ 
dp 



71 

00 du 

■oo 



e~ u2 e u " 



Then 



2% 



2% 



2% 



dp <p {p) du e- u2 e-^ {p) (l - d/dp)e u2d/9p \e* Q <p (p) 
dp <p (p) du e _1l V* 0(p) (l - d/dp) (e^ Q{p+u2) ^ Q {p + u 
dp o (p) du e -^ 2 e \{Q(v+u*)-Q{p)) x 



x 



<p (p + u 2 ) 1 - -Q'(p + u 2 ) - #,(p + it ) 



2i 



dp <p {p) du e ~ u2 e^ Q{p+u2) - Q{p)) x 



x 



— (p + m 2 ) 2 ^ (p + m 2 ) - #>(p + « 2 ) 



Changing variables according to u 2 = Vy/a and p = y\/a, this becomes 



(z) 



2ia± 



dy <p (y) 



dv 



e -v^ e \{Q(y+v)~Q(y)) 



X 



X 



1 _ 1 

-(y + vf<f (y + v)--= <f' (y + v) 
1 Wat 



where Q(y) = y/a (2y - y 3 /3). 

By numerical integration we obtain 



(z) = -i x 0.869 



(41) 



(42) 



(43) 



(44) 



compared with the Schwinger- Dyson estimate — i x 0.908 of Eq. (I2ip . We should also note 
that, when scaled to the value g = 1, our results agree with those of Ref. 15j, apart from 
the fact that the sign of G% should be positive, not negative as it is there. 
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V. DISCUSSION 

We have calculated the Green functions for the wrong-sign quartic by exploiting the 
similarity transformation to an equivalent Hermitian Hamiltonian, and used these as a check 
of the method of truncation of the Schwinger-Dyson equations. We explicitly discussed the 
case m = 0, but the extension to m 7^ is straightforward, as is the extension to higher-order 
Green functions once the identity (j4"Tl) is taken into account. We have also explained why 
previous calculations of these Green functions, which did not explicitly involve the metric, 
nonetheless incorporate it correctly. 

In quantum mechanics it happens that we are able to calculate the metric exactly, but 
the method unfortunately does not extend to higher dimensions. However, for the purposes 
of evaluating Green functions it is not necessary because one can use non-perturbative 
techniques applied to the Schwinger-Dyson equations or the associated functional integrals, 
in the knowledge that such methods are guaranteed to evaluate the correct Green functions. 
As far as the functional integral is concerned, the strong-coupling expansion used in Ref. [I5J 
is a viable method only for the m = case. Unfortunately, Monte-Carlo simulation is not 
viable, because of the complex nature of the functional integral. Altogether it seems that 
the Schwinger-Dyson approach is the most promising one. 
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